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Slow flows of an ideal compressible fluid (gas) in the gravity fleld in the presence of two isen- 
tropic layers are considered, with a small difference of speciflc entropy between them. Assuming 
irrotational flows in each layer [that is Vi^2 = ^^1.2], and neglecting acoustic degrees of freedom by 
means of the conditions div(p(2:) V<^i,2) ~ 0, where p{z) is a mean equilibrium density, we derive 
equations of motion for the interface in terms of the boundary shape z = ri(x, y, t) and the differ- 
ence of the two boundary values of the velocity potentials: '4^{x,y,t) — tlii — ■02. A Hamiltonian 
structure of the obtained equations is proved, which is determined by the Lagrangian of the form 
C — J p{ri)rit4> dxdy — Hlrj, ip}- The idealized system under consideration is the most simple theo- 
retical model for studying internal waves in a sharply stratifled atmosphere, where the decrease of 
equilibrium gas density with the altitude due to compressibility is essentially taken into account. 
For planar flows, a generalization is made to the case when in each layer there is a constant potential 
vorticity. Investigated in more details is the system with a model density proflle p{z) oc exp(— 2qz;), 
for which the Hamiltonian Wjry, tp} can be expressed explicitly. A long-wave regime is considered, 
and an approximate weakly nonlinear equation of the form Ut+auu^ — b[—dx + a^]^^^Ux = (known 
as Smith's equation) is derived for evolution of a unidirectional wave. 

PACS numbers: 47.10.Df, 47.55.-t, 92.60.-e 



I. INTRODUCTION 

Internal waves constitute an important part in the dy- 
namics of such complex systems as are Atmosphere and 
Ocean (see, e.g., Refs. [ll-[lO|. and references therein). 
These waves are known to propagate at the background 
of some inhomogeneity of internal properties of the fluid 
(gas). In the ocean the main role is played by salt con- 
centration and temperature, while in the atmosphere the 
most important factors are specific entropy and air mois- 
ture. Non-uniformity of shear flows should be mentioned 
as well. The internal wave dynamics depends essentially 
on the condition if the stratification is smooth enough in a 
wide range of altitudes, or the change of internal proper- 
ties takes place sharply near some surface. The last case, 
as a rule, is more convenient for a theoretical study, since 
the spatial dimensionality of the problem is reduced. In 
many works therefore simplified atmospheric and oceanic 
models are considered, where the system consists of sev- 
eral layers, with homogeneous fluid within each layer, and 
then the dynamics of interfaces between the layers is in- 
vestigated (see, e.g., [lll - [T7| . and references therein). To 
the best author's knowledge, in all previous finite-layer 
models the fluid was assumed to be incompressible, even 
when the atmosphere was modeled. In the present work, 
perhaps for the first time, an essentially compressible 
two-layer atmospheric model is considered. Here it is as- 
sumed that there is a sharp boundary z = r](x, y, t) sepa- 
rating two regions of potential flow, with a constant value 
of specific entropy in each layer. The relative diff'erence 
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of that values is small, and it ensures the slowness of typ- 
ical flow velocities compared to the local speed of sound. 
Accordingly, the acoustic degrees of freedom can be ef- 
fectively "filtered" by the conditions V • (p(z)v) = in 
each layer (where p{z) is the equilibrium density), instead 
dealing with the full continuity equation pt-l-V- (pv) = 0. 
This idea to eliminate relatively fast sound waves was 
used previously to obtain simplified equations describing 
convection and internal waves in a continuously stratified 
compressible fiuid [18- 20], and also slow isentropic vortex 
fiows in a compressible fiuid placed in a static external 
field [U [ill . The distinction of the present model is that 
the potentiality condition in each layer, together with the 
equation V • {p{z)Vif) = for the velocity potential, al- 
low us to represent equations of motion in terms of the 
interface shape z = 7]{x,y,t) itself and the difference of 
the two boundary values of the velocity potential. More- 
over, we succeeded in proving a Hamiltonian structure of 
the obtained equations, which is a generalization of the 
canonical structure discovered by V. E. Zakharov in the 
dynamics of waves at the free surface of an ideal incom- 
pressible fluid j23l - l25j . In the two-dimensional (2D) case, 
it is possible to consider in the framework of the two- layer 
model also shear flows with piecewise constant potential 
vorticity. The Hamiltonian theory is naturally modified 
in that case. As applications of the developed theory, 
we obtained the dispersion relation for internal waves 
in the two-layer compressible atmosphere, and we de- 
rived a nonlinear equation which is intermediate between 
the Korteweg-de Vries and the Benjamin-Ono equations 
[26I [27| . This equation determines slow evolution of a 
unidirectional wave and it takes into account the dis- 
persive correction of a special form, taking place in the 
model. Previously, a similar equation was derived in a 
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different physical context by Ronald Smith [28|, who in- 
vestigated continental-shelf waves in ocean. 

The paper is organized in the following way. In Section 
2, simplified equations for the two-layer compressible at- 
mospheric model are suggested, and their Hamiltonian 
structure is proved. In Section 3, calculations are per- 
formed for the case of exponential profile of equilibrium 
density, including derivation of the dispersion relation 
and the approximate nonlinear equation for evolution of 
propagating wave. In Section 4, the generalization of the 
model to 2D flows with a piecewise constant potential 
vorticity is made. In Section 5, some conditions of ap- 
plicability of the model are briefly discussed, as well as 
perspectives of future research. Finally, in the Appendix 
we derive three-dimensional (3D) Green's function deter- 
mining the Hamiltonian of the system in the presence of 
the lower flat boundary. 



II. APPROXIMATE EQUATIONS AND THEIR 
HAMILTONIAN STRUCTURE 



Let us assume that in the equilibrium state the first 
layer of gas occupies the region < z < h and has 
the density pi(-z), while the second layer occupies the 
region z > h and has the density P2{z) [for simplicity, 
we have supposed that the lower rigid boundary — "the 
Earth surface" — is fiat, but the more general case of 
nontrivial topography can be considered in analogous 
way]. Of course, functions pi{z) and p2(-z) cannot be 
arbitrary, since in fact they are specified by the hydro- 
static balance condition together with an equation of 
state of the gas (see below). For the further derivation 
of approximate equations describing potential flows in 
this system which are slow compared with a local speed 
of sound c, the following condition is very important: 
(pi-p2)<p(z) = (pi+p2)/2. 

The starting-point equations for potential isentropic 
gas flow in each layer are the non-stationary Bernoulli 
equation and the continuity equation. 



'w{p) — gz + const, 



dtp + V ■ (pVip) = 0, 



(1) 
(2) 



where (p{r, t) is the potential for the velocity field v, satis- 
fying the condition of zero normal derivative at the rigid 
boundary, that is dzip{x,y,0) = 0; p{r,t) is the density, 
w{p) is the specific enthalpy which is defined by the for- 
mula 



wip) = wi,2(p) 



" dpiAp) 



P 



(3) 



Here p = pi^2(p) is the pressure as a function of density 
in each layer, with P2{p) - Pi{p) < [P2{p) +Pi{p)]/'2.- In 
the equilibrium state the velocity potential = 0, the 
enthalpy 1^1,2(^1,2(2;)) = consti^2 — gz, and the pressure 



is related to the density by the hydrostatic formula 
Pi,2{z)^Po-g pi,2{z)dz. 

Jh 



(4) 



Let us consider slow flows when pi,2 = ^1,2(2) +Pi,2 and 
u'1,2 ~ consti^2— gz + pi,2/p{z), where pi,2 are relatively 
small corrections to the pressure field due to fiuid fiow. 
The equations of slow motion in the main order in v/c 
take the form 



(V(^i 



Pia_ 

2 + -p{z) 

V • (p(z)V(^i,2) 



- 0, 



= 0. 



(5) 
(6) 



It is the neglect of time derivative dtp in the continuity 
equation that allows us to exclude from the considera- 
tion acoustic degrees of freedom and retain only "soft" 
modes as the internal waves which are conditioned by the 
relatively small difference of the two equilibrium density 
profiles. Compressibility of the medium in this model is 
manifested in form that a volume of each fluid element 
at slow motion is effectively "adapted" to the equilibrium 
density p{z), expanding when going up and compressing 
when going down [since p'{z) < 0]. 

Let the shape of disturbed interface be given by equa- 
tion z = ?7(x, i), where x = {x,y) is the radius-vector in 
the horizontal plane, and let the boundary values of the 
velocity potentials be ^i,2(x, t) = iy9i.2(x, r]{x, y, t),t). At 
the free interface, the normal component Vn of the veloc- 
ity field should be continuous, as well as the pressure. It 
is also clear that a local speed of boundary motion in the 
normal direction [for definiteness, the normal vector n is 
directed from the first layer to the second one] is equal to 
Vn- From these considerations, two kinematic conditions 
and one dynamic condition are derived, which determine 
evolution of the system: 



d'-P2 
dn 



dn 

m - Vn^T+WW 



f[(Vv'i)^ 



-g I [pi{z) 

Ih 



P2{z)]dz = 0. 



(7) 
(8) 

(9) 



It follows form Eq.© that i/'i and ip2 are related to 
each other by a linear integral dependence. Therefore, 
if wc fix the difference ip(ii.,t) = ipi — ■02, then each po- 
tential will be fully determined. Taking into account the 
equalities 



dipi,2/dt = [d(pi^2/dt + {d(pi^2/dz)rit] 



(10) 



it is easy to check that the equations of motion for the 
two main functions 77(x, t) and ip{x, t) possess the Hamil- 
tonian structure 



piv)vt 



sn 

dip' 



sn 

Srj ' 



(11) 
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with the corresponding Lagrangian 



(12) 



The Hamiltonian functional 'H{?7, V'} is given by the fol- 
lowing expression: 



U 



i)(x) 



dz 







+ 



where 



(13) 



(14) 



that is the Hamiltonian Ji is the sum of the kinetic energy 
and an effective potential energy. Let us prove the above 
statements. 

Indeed, the variation Sip entails some variations Sipi_2, 
and consequently — a variation of the kinetic energy. 
The corresponding variation of the Hamiltonian after in- 
tegration by parts is determined by a surface integral 
along the interface z = ^^(x), and it takes the form 



sn 



Sip 



= / p(Vv3i • n)StpidS - / p(V(^2 • n)Sip2dS 
Js Js 

= I pivWnVl+WWSH^x. (15) 



From here we have ST-L/Sifj — p{rj)Vn^y^ + (Wr|)^ and, 
making comparison with Eq.®, we prove the first equa- 
tion from Eqs. dTT]) . Calculation of variational derivative 
5T-L/5rj is slightly more complicated, because when the 
integration domain is varied, we have to ensure that af- 
ter the interface variation the difference ipi — "02 takes 
at the new boundary the same value V'(x) which was 
before the variation at the old boundary. It is easy 
to understand that due to the above requirement the 
values of the potentials at the place of the old bound- 
ary are changed after variation 5rj by small quantities 
5'tl>f2 — ~^'n{dzfi,2)\z=i)- Accordingly, variation of the 
kinetic energy in this case consists of two contributions. 
The first contribution comes from the change of integra- 
tion domain: 



5rj 



(V^2)' 



Srjd^s.. 



(16) 



The second contribution is related to the changes of the 
potentials ^pi^2 in non-varied domains due to variations 
of their boundary values by the quantities Sipi''!^- It is 
easy to understand that this contribution is equal to 
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dr]d'^x. (17) 



pVn\/l + (V77)2[(9^(^2 -d^ipi] 



Taking into account also variation of the effective poten- 
tial energy, we obtain as the result 



sn 



(V^2)' 



9 



P2{z)]dz 
(18) 



Looking at Eqs.® and pO)) . we obtain from here the 
second equation of the Eqs.(fTT|). 

The Hamiltonian nature of the system under consid- 
eration in principle allows us to apply to it the stan- 
dard set of methods [25[. However, a technical difficulty 
is that the kinetic energy is not expressed directly but 
through solutions of the partial derivative equation 
with non-constant coefficients, and in domains with a 
curved boundary z = ry(x). Let us nevertheless suppose 
that particular solutions of Eq.([B]) are known in the form 
of linear combinations 



(pk(x,z) = [A<i>(-)(z) 



B^i+\z)y 



(19) 



with decaying at z — > +oo functions $^ (z), and with 

growing at z — )■ +oo functions ^\^\z). In other words, 
for every k a general solution is known for the following 
equation, 



$"(z) + ^$'(z)-fc2$(z) = 0. 



(20) 



Then for approximate calculation of the Hamiltonian at 
small deviations C(x, i) — '7(x, t) — h, with the condition 
|VC| ^ 1, one can write 



(pi(x, z) 
'y32(x,z) 



(2^)2 



l(2)^(-) 



(27r) 



2<^n 



iz)e 



'(z)]e'''-(21) 
(22) 



After that from the set of boundary conditions 1) 
9,(^1 (x,0) = 0; 2) 5„(^i(x,/i + C(x)) = 9„<^2(x, /i-f C(x)); 
and 3) ipi{x,h + C(x)) - '^2(x, /i -I- C(x)) = '(/'(x) it is 



(1) 



and 



possible to find the unknown function Aj^ 

(2) 

Ajj. [and consequently the required quantity y„(x)] in 
the form of an expansion in C. Such a method of present- 
ing the Hamiltonian as a series in the small parameter 
of characteristic wave steepness is generally used in the 
theory of surface water waves |23l425| . In particular, this 
method allows us to obtain the dispersion relation for 
low-amplitude internal waves: 



Di{h,k)D2{h,k) 
[D2{h,k)+Di{h,k)y 



(23) 



where g{h) is a renormalized gravity acceleration: g{h) ~ 
g[pi{h) — p2{h)]/p{h), and the short-hand notations have 
been used: 



Di{h,k) 



'W(/.)$^)(0) 



,(24) 
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(25) 



Note that Di(/i, fc) > and D2{h, k) > 0. 

As to the system under consideration, here in some 
cases another way can be suitable how to calculate the 
Hamiltonian. Since the kinetic energy takes the form 



= ^ y (j • ^)d^^dz, 



where j = pv is the divergence-free field of the current 
density, we can introduce for j a vector potential A which 
satisfies the equation 



curl _ curl A. = ^l — curl v. 

p(z) 



(26) 



with the boundary condition [dxA'^y\x,y,Qi) — 
dyA^^^{x,y,G)] = 0. After that the kinetic energy 
can be re- written as follows, 



/C = -jA-n cPxdz 

= iy"G,fc(ri,r2)f],(ri)nfc(r2) d^nd^ra, (27) 

where Gik{ri, r2) is the Green's function for Eq. ((26|) . As 
far as the (singular) vorticity field ft is totally concen- 
trated at the interface z = 7?(x), and the vortex lines 
coincide with levels of the function tp(x.) at that surface, 
the half-space integration will reduce to integration along 
the surface z — rjix) by means of the change 

(28) 

As the simplest example, in this work an exponen- 
tial profile p{z) = po exp(— 2az) of the equilibrium den- 
sity will be considered, when Ea. (PE| after substitution 
A = poe~^"^F turns into an equation with constant coef- 
ficients. Generally speaking, if taken globally, such a de- 
pendence contradicts to adiabatic equations of state for 
real gases, for those we rather have p Cip'^ , where 7 
is the adiabatic exponent [for single-atom gases 7 = 5/3, 
for gases consisting of two-atom molecules 7 — 7/5], and 
therefore p(z) « C2(zo — zy^'''^~^\ where zq is the alti- 
tude of the upper edge of the atmosphere. Nevertheless, 
locally on the vertical coordinate near z = h, every re- 
alistic dependence p{z) is approximated by an exponent, 
provided not very long waves are considered. We still 
would like to note that the case p(z) « C2{zq — z)-'^/^'''"^^ 
also admits analytic investigation, though more difficult, 
since the functions $[.^''(z) in the particular solutions 
of Eq.® are expressed in that case through the 
modified Bessel functions and K,,, with the index 



i/- [(7-1)- 



l]/2: 



$1 \z) = [kizo-zT'^I.ikizo-z)), (29) 
^i+\z) = [k{zo-z)]-^K,{k{zo-z)). (30) 



III. THE CASE OF EXPONENTIAL PROFILE 
OF EQUILIBRIUM DENSITY 

Thus, we have to find the Hamiltonian of our system 
in an explicit form for p{z) = poexjp(—2az), and at the 
beginning we will solve Eq. ((26| . Consider here simpler 
case ah 3> 1, when the presence of the flat lower bound- 
ary at z = is not important, because the correspond- 
ing contribution will be shown later to be of the order 
exp(— 2/ia). More cumbersome 3D solution for the vec- 
tor potential in the presence of the boundary z = is 
given in the Appendix. To solve Ea. (PB|) . we use the 
substitution A = poe~'^"^F and re-write the equation in 
Fourier representation: ik x [(zk — 2aez) x Fk] — Ok- 
Applying the well-known formula for the double vector 
cross-product and choosing the gauge (k • Fk) — 0, we 
immediately arrive at a simple equation 



[fc2 + 2ia(k ■ e^)]Fk = rJk- 



(31) 



Now we write down the decaying at the infinity solution 
of the above equation: 



F(r) = 



d^k 



(27r)3 [fc2 + 2ia{k ■ e,)] 
exp[a(z — zi — |r — ri| 



47r|r — ri | 



-n{ri)d^ri. (32) 



Accordingly, the kinetic energy of the 3D system, without 
taking into account the fiat rigid boundary, is given by 
the following expression: 



p-Q|r2-ri I 



g-a(22+^i)f2(r2) • n(ri)dhid^r2. 

\r2 - ril 

(33) 

Passing with the help of formula (|28l) from the space inte- 
gration to the surface integration where singular vorticity 
field is distributed, we arrive at the expression in terms 
of r] and ip, 



^ ^ Po /"exp[-av/|xi-X2p + {qi-mY - a(?7i+'72)] 



87^7 ^|xi-X2|2 + (771-772)2 

X{V7^1 • V1P2 + [V7/>i X Vrii] ■ [V7/>2 X V7y2]} 

xd^x^d^x^, (34) 

where V77 and Vt/j are 2D gradients. If necessary, a 
weakly nonlinear regime in the wave dynamics can be 
easy considered through expansion of the above expres- 
sion in powers of and C,. 

Let us now turn our attention to planar flows. Note 
that in 2D case Fk and ilk are in the essence (pseudo) 
scalar quantities. The presence of the boundary at z = 
can be taken into account by a variant of the "image 
method" , and as the result we have 

1 



F{x,z) = — 
zvr 



ifo (a ^ [x - xiY + {z - zif) 

^Ko{a^{x-XiY + [z + ZiY) 
xe"(^-^i'l7(a;i,zi)da;idzi, (35) 
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where Ko{r) is the well-known Macdonald function. We 
provide below two of many possible integral representa- 
tions for this function: 



^2j^ g'ifeia+ifc2 6 



2n kf + ki + 1 



+00 



exp(ifca — + 1) 



1 



dk. 



(36) 



Consequently, the Green's function in this case takes the 
form 



[ct^J {xi - X2Y + {zi + Z2Y^ 



-0(21+22) 



(37) 



The expression for the kinetic energy of the two-layer 
flow looks as follows: 



/C 



2D 



P0_ 

in 



[xi - X2f + {r]x - 772)^) 



-Kq [a^ (xi - X2Y + (?7i + n-ip- 



(38) 



where V'' — dijj/dx. Use of formulas (l36l) allows us to 
represent this functional in a slightly different form: 



the small arguments, we obtain an approximate kinetic 
energy functional up to the first order in ar], 

^*{'nM = yJ Va-2av){i:'fdx 

{ib'rj)[-dl+a^]^^\ij'r])dx. (42) 



Let us introduce a new unknown function q = [1 — 
exp(— 2a77)]/(2Q!), which up to the constant factor pQ is 
the canonically conjugate for function tjj, and then re- 
write the approximate Hamiltonian in terms of q and ip: 



^ I q{l-aq){^P'fdx 



Pa 



2 

Poff(O) 



(VA)h(9^ + a2]i/2(^'<7)da; 

3" 



2 



a/3 



q 



dx, 



(43) 



where /3 is a dimensionless parameter depending on be- 
havior of the difference [pi{z) — P2{z)] near z = 0. Con- 
sidering propagation of relatively small but finite dis- 
turbances q{x,t) = q{x,t) — q, it is possible by a stan- 
dard procedure to derive weakly nonlinear equation for 
u{x,t) — ipx, which describes a slow evolution of uni- 
directional wave under the influence of weak dispersion: 



2D 



PO ^ 

Q-\Vl-V2\Vk'^+a'-' 



-(vi+V2)Vk'-'+a'-'-\ 



2Vk^ 



2n 



(39) 



As it will be shown later, such a representation is suit- 
able for consideration of long-wave asymptotics in the 
nonlinear wave dynamics. Besides that, it also allows 
us to find easily the dispersion relation for linear waves. 
Indeed, from Eq. (j39p it is obvious that in the quadratic 
approximation the Hamiltonian is given by the formula 



U 



[2] 
2D 



Poe' 



-2ah 



\l 



2Vk^+^ 



-k ip-ki'k 



+9{h)C-kCk 



dk 
2^' 



(40) 



Solving the corresponding linearized equations of motion 
for the Fourier components Cfc(^) ''Pk{t), we find quite 
nontrivial expression for the dispersion relation: 



-2hVk'^+a'^] 



2VF 



(41) 



Note, the same dispersion law takes place in the 3D case, 
due to the isotropy of the system in the horizontal plane 
[it is also confirmed by the formula ([231)] • 

Now we consider the limiting case arj <^ 1 and typical 
wave numbers k satisfying the conditions arj ^ fcr/ <C 1. 
Expanding the exponents in integral (j39l) in powers of 



Ut + cux + auux 



v2ll/2 _ 



= 0, (44) 



where the speed of long linear waves is c w [5(0)9]^^^, and 
the coefficient d w 3/2. Equation of such kind is called 
sometimes "Smith's equation" after the work by Ronald 
Smith [2^] where it arose for the first time in context of 
continental-shelf oceanic waves. It is interesting to note 
that the special form of the dispersive term makes the 
above equation intermediate between the two famous in- 
tegrable models, namely the Kortewe g-de Vries equation 
and the Benjamin-Ono equation 2^ 27 1. In this sense 
the Smith's equation is similar to the Intermediate Long 
Wave equation (ILW) (see, e.g., [IMl [H, H^l ) , but con- 
trary to ILW the Smith's equation is not integrable, as it 
was estabhshed in Ref . [3111 . 



IV. PLANAR FLOWS WITH PIECEWISE 
CONSTANT POTENTIAL VORTICITY 

Now we would like to make an important generaliza- 
tion of the Hamiltonian theory which is possible for 2D 
isentropic flows [in (x, z) plane], namely we will take into 
account the fact that potential vorticity 7 = — fj'^'/p in 
the 2D case is governed by the advection equation 



7t + V • V7 = 0. 



(45) 



This conservation law for the potential vorticity along 
each fluid particle trajectory allows us at consideration of 
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planar flows with a piecewise constant function j^x, z,t) 
to follow only the motion of boundaries where 7 is discon- 
tinuous. In the present paper it is assumed for simplicity 
that 7 has a single jump, and this jump coincides with 
the interface between the layers z — ri[x,t), but gener- 
ally this coincidence is not necessary and a separate curve 
z = ri^,{x,z,t) can be considered where the jump takes 
place. 

Let (sufficiently small) potential vorticities in the lay- 
ers be 71, 2j so that the corresponding stationary shear 
flows C/i,2(z) <^ c satisfy the conditions (we neglect the 
difference between pi_2 and p) 



(46) 



We shall suppose that in the stationary state the veloc- 
ity profile has a "break" a.t z — h, that is C/i,2(-z) = 
-7i,2Ai(z), where 



m(2) = / (47) 

Jh 

A 2D velocity field in each layer now takes the form 

vi,2(a;,2:,t) = {Ui^2{z) + d.j;ipi^2{x,z,t), dzipi,2{x, z,t)), 

(48) 

with the potentials ipi^2 satisfying the same equation ([6|): 
V • pV(/?i,2 — 0, and it implies the existence of the corre- 
sponding stream functions 'i?i,2(a;, z,t): 



pdx^Pia = dz"di.2, pdzipi.2 = -dx-ffi 



(49) 



Instead of Eq. ([5]) , we have to deal now with its general- 
ization: 



dtVl,2 + 7l,2l9l,2 + Ui,2{z)dx(pi,2 + 



(Vyi.2)^ Pi. 2 

2 p{z) 



(50) 

which regards the 2D Euler equation in the case of con- 
stant potential vorticity under the condition V • (pv) = 0. 
Taking into account that the full stream functions of the 
flows under consideration are 

61,2(2;, z, t) = i?i,2(a;, z, t) - t/2 2(z)/(27i,2), (51) 

equation ([50)) can be also represented as follows, 



Ot<yfl,2 + 71, 261,2 H h — -- = 0. 

2 p{z) 



(52) 



Now we note that at the interface z — ri[x,t) there are 
the equalities 



-da:'di{x,'q{x)) = ~dxQ2{x,vix)) ^ p{v)vt 



(53) 



where Vn — {"vi ■ n) ~ (v2 • n). 

Demanding the pressure field to be continuous at z = 
r]{x, t) and reasoning analogously to the case 71,2 = 0, we 



conclude that the evolution equations for the 2D system 
possess the following structure, 

p(77)?7t = 5n/5^, (54) 
-p{r])^l}t+ip{'n)d^^[p{ri)-qt] = SU/Sri, (55) 

where 7 = (71 — 72), and the Hamiltonian % is equal to 
the sum of total kinetic energy and the effective potential 
energy. By a direct calculation it is easy to check that 
the corresponding Lagrangian for the above equations is 

£ = / ^ljptdx+'^ pd^^ptdx ~'H{p,'ip}, (56) 



where p = p{j]) [see Ea.(|47p]. For internal waves in an in- 
compressible liquid, an analogous structure was obtained 
in Refs.fTsl. [l7j. with the difference that in our case p{j]) 
is a nonlinear function (see also Ref.js^] about waves at 
the free surface of a 2D incompressible fluid with a con- 
stant vorticity). 

It is interesting to note that in the quadratic approxi- 
mation the Lagrangian (j56p take the form 



£[2] = -p{h) / ^Ct dx + 



ip\h) 



(57) 

Moreover, it is easy to show that the functional H^^' {C, "0} 
does not depend on 71 and 72 [dependence on 71 and 72 
appears only in higher orders] : 



-p{K) 



dk 

[Nih,k)e^p_k^k + ~gih)C-kCk] ^■ 



(58) 

Function N{h, k) is expressed through the Green's func- 
tion G[{x2 — a;i), zi, Z2] by the following formula: 



j5{h)N{h,k) = 



+00 



G[x, h, h]t 



-'''''dx. 



(59) 



It should be noted that ujQik) — g{h)k^N{h, k) is the dis- 
persion law in the case 71 = 72 = [compare with (|23p ]. 
For example, with the exponential proflle of the equi- 
librium density the quadratic Hamiltonian is given by 
expression (j40p . Solving the corresponding linear equa- 
tions. 



Ck = N{h,k)k^iJ;k, 



-■>Pk + ip{h)r 



Ck 



k 



(60) 
(61) 



we obtain the dispersion law for linear waves at 7 ^ 0: 
1 



-jp{h)kN{h,k) 



+ \j Ylp{h)kN{h, k)f /4 -I- g{h)k^N{h, k). (62) 

Since the singular part of the vorticity field (concen- 
trated at the interface) is determined by the relation 
— fis = W — lp{fi)\^{z — ?7(a;,i)), where 5{z — rj{x,t)) 
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is the Dirac's function, it is convenient to introduce the 
new unknown variable, 

= ?A - 7a~V- (63) 

In variables {fJ,,p} the Lagrangian takes the form (the 
sign in front of the second term has changed) 

£ = Jpfitdx-^ J fid^^^fitdx - 'H{tJ,,p}. (64) 

Now, besides the singular part of the vorticity, there is 
also a distributed part, and the full vorticity field is given 
by the formula 

-n{x,z,t) = p'{x,t)S[z - r]{x,t)] + 

+ l2Piz)+lp{z)e[i]ix,t)- z], (65) 

where 9[r](x,t) — z] is the unit jump function (the Heav- 
iside's function). The Hamiltonian of the 2D system 
is determined with the help of the Green's function 
G[{x2 — xi), zi, Z2] by the following expression: 

n ^ gj Wirj)dx + ^ J G[ix2~xi),zi,Z2] 

xfl(xi, zi)n{x2, Z2)dxidzidx2dz2, (66) 

where Eq. (p5|) should be substituted, and after the in- 
tegrations 77 should be expressed through fi. Let us re- 
mind that in the case p{z) = po exp(— 2a2;) the Green's 
function is given by Eg . ([BT)) . Let us also note that 
in the absence of the density jump a class of flows is 
possible with p = 0. In that case the dynamics of 
the vorticity waves is determined by the Lagrangian 
C-y = —(7/2) J pd~^ptdx — H-ylp}, and the dispersion 
law for such waves is expressed by the formula (j62p where 
g{h) = should be put. 

V. DISCUSSION 



that in the long- wave limit (in the Earth conditions it cor- 
responds to hundreds and thousands kilometers), nonuni- 
form horizontal motions of the whole atmosphere become 
important. Those flows are approximately described by 
a "shallow water theory" with adding the Coriolis force, 
and they lead to variations of a quasi-equilibrium den- 
sity profile. Besides that, the Coriolis force violates the 
potentiality of the flow. Thus, the suggested here theory 
can describe waves with lengths not longer than a few 
kilometers 

In the present work, only first steps have been made in 
the study of internal waves in the atmosphere within the 
compressible two-layer model. Promising directions of 
further research can be outlined as follows. First, a gen- 
eralization of the model is evident for more layers and for 
continuous limit, which will enrich it because an interac- 
tion between several interfaces in many cases is able to 
introduce new interesting effects as instabilities etc. Sec- 
ond, we should mention a wide class of problems about 
interaction of internal waves and mountains, which also 
can be studied with the help of this model. Third, nonlin- 
ear wave dynamics can be simulated numerically. Fourth, 
an analogous Hamiltonian formulation is possible for con- 
sideration of axisymmetric flows with a piecewise con- 
stant generalized potential vorticity. Fifth, it seems likely 
that analogous finite-layer models are possible not only 
in the Eulerian hydrodynamics, but in a wider class of 
conservative hydrodynamic systems as well, for instance, 
in the hydrodynamics of a relativistic fiuid placed in a 
strong static gravitational field described by a metric 4- 
tensor. Accordingly, there is a perspective of application 
of a similar theory to astrophysical problems, where the 
equilibrium density possesses the spherical symmetry, as 
a rule. 

These investigations were supported by RFBR (grants 
09-01-00631 and 07-01-92165), by the "Leading Scientific 
Schools of Russia" grant 6885.2010.2, and by the Pro- 
gram "Fundamental Problems of Nonlinear Dynamics" 
from the RAS Presidium. 



In this work, a compressible two-layer atmospheric 
model has been suggested, intended for theoretical study 
of internal waves at the interface between two isentropic 
layers of a gas with nearly equal values of specific en- 
tropy. In the derivation of the approximate equations 
it was supposed that the fiow velocities are small com- 
pared with a local speed of sound. It should be noted 
that this condition puts the lower limit for characteris- 
tic wave numbers: k ^ |p'(/i)|/p(/i), because at longer 
scales the velocity field penetrates into the upper layer 
rather far, where in view of constant entropy the temper- 
ature is small together with a local speed of sound, and it 
violates the starting-point assumption of the model. To 
some extent the above limitation is softened if somewhere 
above the second layer there is the third layer, with very 
high temperature, and therefore the boundary between 
the second layer and the third layer can be effectively 
treated as a "rigid lid" . However one should remember 



Appendix A: Correction to 3D Green's function due 
to the flat boundary 

To satisfy the boundary condition [dxF^y'>{x,y,0) — 
dyF'-'=^x,y,0)] = 0, which ensures zero normal com- 
ponent of the velocity field at the rigid flat bound- 
ary, we add to the particular solution (15^ of the non- 
homogeneous equation (j31[) some specially selected solu- 
tion of the corresponding homogeneous equation, decay- 
ing at z — > +00: 

F(-)(x,z)= / 4%fkexp[ik-x + z(a- v/fc^T^)], 
J (27r)^ 

(Al) 

where fk ~ ifk\ fu \^) satisfles the condition of 2D 
transversal gauge (k • fk) = 0. It is not difficult to un- 
derstand that fk should be taken in the following form 
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(here and later on k and v are tensorial indices in the 
horizontal plane): 



27r (fc2 + ^2 ^ 2ia^) 



(A2) 



where f7('')(k,^) = / ^^(''^(xi, zi)e-*-^i-^«^id2-^^^^^ jg 
the Fourier image of the horizontal component of the 
vorticity field. Now we transform the integral (jA2l) : 



r!H(xi,zi)e 2 

2^ (fc2 + e2 + 2mO ' ' 



2V^ 



(i2xidzi(A3) 



For ^-integration we have used the fact that r2(x, z) 
is non-zero only at z > 0, and therefore the integra- 
tion contour can be closed in in the lower complex half- 
plane. Then we substitute the result into Eq. ljAip and 
see that the conditioned by the flat boundary correction 
Gii/^(xi, X2, zi, Z2) to the Green' function actually de- 
pends on the variables x = X2 — Xi and s = zi + Z2 , and 
it is expressed by the following formula: 



G'L^(x,s) = Pot 



fc2 



2\/fc2 



(27r) 



~as,^ o o X-1^6Xp[-a\/x2 + s2] 



where ^ is the inverse 2D Laplace opera- 
tor. Let us introduce the notation Z?(|x|,s) = 
A^-'- [exp(— Q!\/x2 -I- s2)/v'x2 -I- s2]_ jjj virtue of the 
definition, function D{r, s) satisfies the equation 



I d ( dD\ _ exp(-aVr2 + s2) 
-rd-ry^)" 7^T^ ' ^^^^ 

from which we obtain by a simple integration 



rdD /dr — [exp(— as) — exp(— r'^ -t- s^)]/a. (A6) 

It should be noted that the second derivatives 
df^di,D(\K\, s) can be expressed through the combination 
r~^dD/dr: 



a.-,9.D(|x|,s) = 



, , x^x, d (ldD\ , ^ IdD 
r or \r or J r or 



(A7) 



Collecting the obtained expressions and taking into ac- 
count Eq. (jA6p . we write the required correction to 3D 
Green's function in the final form: 



G(;)(x,s)-po is.. -2 



Xj^Xi/ 



[exp(— 2q;s) — exp(— as — a\/yp- -|- s2)] 



47rax2 



+Po — 
V X 



Xv ^ \ exp(-as - a\Jy? -V s^) 

'Oku] — , „ , „ • (A»j 

^ 47rVx2 + s2 
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